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Abstract 

We report the first observation of radiative decay B° — > (j)K 7 using a data sample of 772 x 10 6 
BB pairs collected at the T(45 r ) resonance with the Belle detector at the KEKB asymmetric-energy 
e + e~ collider. We observe a signal of 35 ± 8 events with a significance of 5.4 standard deviations 
including systematic uncertainties. The measured branching fraction is B(B° — > (j)K°^f) = (2.66 ± 
0.60 ± 0.32) x 10~ 6 . We also precisely measure B(B+ (f)K + j) = (2.34 ± 0.29 ± 0.23) x 10~ 6 . 
The uncertainties are statistical and systematic, respectively. The observed M^k mass spectrum 
differs significantly from that expected in a three-body phase-space decay. 

PACS numbers: 14.40.Nd, 13.25.Hw, 11.30.Er 
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Rare radiative decays of B mesons play an important role in the search for physics 
beyond the standard model (SM) of electroweak interactions. These flavor changing neutral 
current decays are forbidden at tree level in the SM, but allowed through electroweak loop 
processes. The loop can be mediated by non-SM particles (for example, charged Higgs or 
SUSY particles), which could affect either the branching fraction or the time-dependent CP 
asymmetry. 

The current measured inclusive world average branching fraction for B — > ((3.55 ± 
0.26) x 10~ 4 (H), is one standard deviation (a) higher than the SM prediction at next-to- 
next-to-leading order (NNLO) (3.15±0.23) x 10~ 4 j2j, and still allows significant new physics 
contributions to radiative B decays. Exclusive b — y 57 decays have also been extensively 
measured, but their sum so far accounts only for 44% of the inclusive rate. Therefore, 
further measurements of branching fractions for exclusive B — y <j)K / -f modes will improve 
our understanding of the b — y 57 process. The neutral mode B° —> <pK 0r y Q can be used 
to study time-dependent CP asymmetry, which is suppressed in the SM by the quark mass 
ratio (2m s /mb) 0, HJ. In several models beyond SM, the photon acquires an appreciable 
right-handed component due to the exchange of a virtual heavy fermion in the loop process, 
resulting in large values of time-dependent CP asymmetries. Due to the narrow width of 
the resonance, the decay B — y (pK'y is well separated from the background and can be 
effectively used for measurements of photon momentum over a wide interval. In addition, 
this mode can also be used to search for a possible contribution from kaonic resonances 
decaying to <f>K. Furthermore, we can probe the photon polarization using the angular 
distributions of the final state hadrons j6|, 0] • 

The decay B° —y (pKg'y can be described by the conventional radiative penguin diagram 
with the creation of an additional ss pair as shown in Fig. [TJ 




FIG. 1: Feynman diagram for the radiative penguin penguin decay B° —y (pK^ with ss pair 
creation. 

The branching fractions for B — y <j)K , ~f decays have already been reported by the Belle 
and BaBar collaborations. Belle measured B(B + -»■ 4>K + ^) = (3.4 ± 0.9 ± 0.4) x 10~ 6 and 
B{B° -> (f)K°j) < 8.3 x 10" 6 at the 90% confidence level (C.L.) using 96 x 10 6 BB pairs @. 
BaBar measured B(B + ->■ <pK + ^) = (3.5±0.6±0.4) x 10~ 6 and B(B° <j>K°i) < 2.7xl0" 6 
at the 90% C.L. using 228 x 10 6 BB pairs jij. BaBar also reported the direct CP asymmetry 
for B^ — y (pK^^f, Acp = (— 26±14±5)%. We report herein the first observation of radiative 
decay B° — y <f)K 0/ ~f and an improved measurement of B + — y <f)K +/ ~f using a data sample of 
772 x 10 6 BB pairs collected at the T(4S) resonance with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider [Toj]. This data sample is nearly eight times larger than the 
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sample used in our previous measurement 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of- flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside the coil is instrumented to detect K° L mesons and to identify muons 



(KLM). The detector is described in detail elsewhere [llj]. Two different inner detector 
configurations were used. For the first sample of 152 x 10 6 BB pairs, a 2.0 cm radius 
beampipe and a 3-layer silicon vertex detector (SVD-I) were used; for the latter 620 x 10 6 
BB pairs, a 1.5 cm radius beampipe, a 4- layer silicon vertex detector (SVD-II), and a small- 
cell inner drift chamber were used. A GEANT-based simulation of the Belle detector is used 
to produce signal Monte Carlo (MC) [l2| event samples. 

The signal is reconstructed in the decays B + — > <j)K +/ -f and B° — > (pKg'-f, with <fi — > K + K~ 
and K s — > tt + tt~. All the charged tracks used in the reconstruction (except for charged 
pions from K^s) are required to satisfy a requirement on the distance of closest approach 
to the interaction point (IP) along the beam direction, \dz\ < 5 cm, and in the transverse 
direction, dr < 2 cm. This eliminates poorly reconstructed tracks or tracks that do not 
come from the interaction region. Charged kaons are identified using a likelihood ratio 
C(K/n) > 0.6, based on information from the ACC, TOF and CDC (dE/dx) detectors. This 
requirement has an efficiency of 90% for kaons with a 8% pion fake rate. A less restrictive 
likelihood ratio requirement C(K/n) > 0.4 is applied to the kaon candidates, which are used 
to reconstruct the meson. The invariant mass of the <fi candidates is required to be within 
—0.01 GeV/c 2 < M K + K - — < +0.01 GeV/c 2 , where denotes the world-average 
mass [l3| . 

Neutral kaon (Kg) candidates are formed from the ir + 7r~ combinations with invariant 
mass in the range 0.482 GeV/c 2 < M w+W - < 0.514 GeV/c 2 . The selected candidates must 
pass a set of momentum-dependent requirements on impact parameter, vertex displacement, 
mismatch in the z direction, and the direction of the pion pair momentum as described in 



the Ref. 14 



The primary signature of this decay is a high energy prompt photon. These are selected 
from isolated ECL clusters within the barrel region (32° < 6 1 < 129°, where Q 1 is the polar 
angle of the photon in the laboratory frame) and center-of-mass system (cms) energy (E™ s ) 
in the range 1.4 to 3.4 GeV. The selected photon candidates are required to be consistent 
with isolated electromagnetic showers, i.e., 95% of the energy in an array of 5 x 5 CsI(Tl) 
crystals should be concentrated in an array of 3 x 3 crystals and should have no charged 
tracks associated with it. We also suppress the background photons from vr°(?7) — > 77 using 
a likelihood C n o(C v ) < 0.25, calculated for each photon pair consisting of the candidate 



photon and any other photon in the event [15 



We combine a meson candidate, a charged or neutral kaon candidate and the radiative 
photon to form a B meson. The B candidates are identified using two kinematic variables: 
the energy difference AE = — £J^™ m and the beam-energy-constrained mass M^ c = 
V (-^beam) 2 ~ (P^Y S where -Eboam is the beam energy in the cms, and E™ s and p c ™ are 
the cms energy and momentum, respectively, of the reconstructed B candidate. In the M^ c 
calculation, the photon momentum is replaced by (E^ m — E^ s ) to improve resolution. The 
events that satisfy the requirements M^ c > 5.2 GeV/c 2 and \AE\ < 0.3 GeV (defined as 
the fit region) are selected for further analysis. Using MC simulations, we find nearly 12% 
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(3%) of events have more than one B candidate for the B + — > cf)K +r ) (B° — > cf)K 0r y) mode. 
In case of multiple candidates, we choose the best candidate based on a series of selection 
criteria, which depend on a x variable using the candidate's mass (and the K° s mass m 
the neutral mode) as well as the highest E^ ms and the highest £(K/ir) in the charged mode. 
For events with multiple candidates, this selection method chooses the correct B candidate 
for the B + (pK + -f (B° (f>K°j) mode 57% (69%) of the time. We define the signal 
region as 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and -0.08 GeV < AE < 0.05 GeV. The AE 
signal region is asymmetric in order to include the tail in the lower region due to photon 
energy leakage in the ECL. 

The dominant background comes from e + e~ — >■ qq (q = u, d, s, or c) continuum events. 
We use two event-shape variables to distinguish the spherically symmetric BB events from 
the jet-like continuum events. A Fisher discriminant [l6[ is formed from 16 modified Fox- 
Wolfram moments [l7[ and the scalar sum of the transverse momenta. The second variable is 
the cosine of the angle between the B flight direction and the beam axis (cos#b) in the cms 
frame. For each variable, we obtain the corresponding signal and background probability 
density functions (PDFs) from large MC samples. A likelihood ratio lZ s /b = C s / (C s + Cb) is 
formed, where C s (Cb) denotes the product of Fisher discriminant and cos 9 b PDFs for the 
signal (background). The selection criteria on TZ s /b are determined by maximizing the figure 
of merit, N$/y/Ns + Nb, where iVg (Nb) is the expected number of signal (continuum) 
events in the signal region. We require TZ s /b > 0.65, which removes 91% of the continuum 
while retaining 76% of the signal. 

In addition to the dominant continuum background, various BB background sources are 
also studied. In the B° — > <pK^ mode, some backgrounds from b — > c decays, such as D°ir , 
D°r] and D~p + peak in the M^ c distribution. We remove the dominant peaking backgrounds 
by applying a veto to (j)Kg combinations consistent with the nominal D mass [Hj]. Some 
of the charmless backgrounds, where the B meson decays to 0fT(892), (pKn and <f)Kr] 
also peak in M^ c . In these charmless modes, one of the photons from a ir° or 77 may not be 
detected in the calorimeter while the other is reconstructed as the signal high-energy photon. 
Therefore, these backgrounds shift towards lower AE. Another significant background is 
non-resonant B — > K + K~ K^y, which peaks in the AE-M^ C signal region. The fraction of 
such events is estimated to be (12.5 ± 6.7)% using the mass sideband, 1.05 GeV/c 2 < 
M K+K - < 1.3 GeV/c 2 , in data. 

The signal yield is obtained from an extended unbinned maximum-likelihood fit to the 
two-dimensional AE-M^ distribution in the fit region. We model the shape for the signal 
component using the product of a Crystal Ball line shape [l8[ for AE and a single Gaussian 
for Mb c . The continuum background is modeled with a product of first order Chebyshev 
polynomial for AE and an ARGUS [HJ function for M^ c . The b — > c background is modeled 
with a product of second order Chebyshev polynomial for AE and an ARGUS plus Gaussian 
function for M^ c . The small charmless backgrounds (except the non-resonant component) 
are modeled with a functional form that is the product of two Gaussians for AE and with 
a single Gaussian for M^ c [2(|. In the final fit the continuum parameters are allowed to 
vary while all other background parameters are fixed to the values from MC simulation. 
The shape of the peaking backgrounds are fixed to that of signal in M^ c and AE. In the 
B + — > <j)K +/ -y mode, the non-resonant background yield is fixed to the value from the <p 
sideband and assuming isospin symmetry, the same non-resonant fraction is used in the 
neutral mode. The signal shapes are adjusted for small differences between MC and data 
using a high statistics B° — > i^*(892)°(— > K + tt^)^ control sample. The invariant mass of 
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the K* candidates are required to satisfy 0.820 GeV/c 2 < M K+7T - < 0.970 GeV/c 2 . The fit 
yields a signal of 136 ±17 B + — > <pK + '~f and 35 ±8 B° — > <pK^ candidates. The projections 
of the fit results onto AE and Mb c are shown in Fig. |2j The signal significance is defined 
as \/—2 ln(£o/£ m ax), where £ max is the maximum likelihood for the best fit and Co is the 
corresponding value with the signal yield fixed to zero. The additive sources of systematic 
uncertainty described below are included in the significance by varying each by its error and 
taking the lowest significance. The signal in the charged mode has a significance of 9.6 a, 
whereas that for the neutral mode is 5.4 a. 
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FIG. 2: The AE and Mt, c projections for B + — > (j)K + ^ (upper) and B° — > </>ifg7 (lower). The 
points with error bars represent the data. The different curves show the total fit function (solid 
red), total background function (long-dashed black), continuum component (dotted blue), the 
b — > c component (dashed-dotted green) and the non-resonant component as well as other charmless 
backgrounds (filled magenta histogram). 

We also examine the (f)K invariant mass distribution of the signal. To unfold the Mal- 
distribution, we subtract all possible backgrounds and correct the <pK invariant mass for 
the efficiency. The background-subtracted and efficiency-corrected M^x distributions are 
shown in Fig. [3J Nearly 72% of the signal events are concentrated in the low-mass region 
(1.5 GeV/c 2 < M^k < 2.0 GeV/c 2 ). It is clear that the observed <j)K mass spectrum differs 
significantly from that expected in a three-body phase-space decay. The MC-determined 
reconstruction efficiencies (defined as the ratio of signal candidates passing all selection 
criteria to the total number of events generated) are corrected for this M^k dependence. 
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TABLE I: The signal yields, significances, weighted efficiencies and branching fractions for the 
B + — > (j)K + j and B° — > (f)K 0/ y decay modes. 



Decay mode 


Yield 


Significance (a) 


Efficiency (%) 


Branching fraction (10 6 ) 


B+ -> <pK+>y 


136 ± 17 


9.6 


15.3 ±0.1 


2.34 ±0.29 ±0.23 


B° -> <pK°~f 


35 ±8 


5.4 


10.0 ±0.1 


2.66 ±0.60 ±0.32 



From the signal yield {N sig ), we calculate the branching fraction (£>) as N sig / (e x N BB - x B sec ), 
where e is the weighted efficiency, N BB - is the number of BB pairs in the data sample, and 
£> sec is the product of daughter branching fractions 
Table E 



13 . The results are summarized in 



500 




1.5 2 2.5 3 3.5 4 

M((|> K) (GeV/c 2 ) 




1.5 2 2.5 3 3.5 4 

M((|> K) (GeV/c 2 ) 



FIG. 3: The background-subtracted and efficiency-corrected <f)K mass distributions for B + — > 
(j)K + ^/ (left) and B° — > <j)Kgj (right). The points with error bars represent the data. The yield in 
each bin is obtained by the fitting procedure described in the text. The three-body phase-space 
model from the MC simulation is shown by the circles (blue) and normalized to the total data 
signal yield. 

We fit the data with each fixed parameter varied by its ±1 a error, and then the quadratic 
sum of all differences from the nominal value is assigned as the systematic error on the 
signal yield. We checked for possible bias in the fitter by doing ensemble tests with MC 
peseudo-experiments. The statistical errors obtained from our measurements are within 
the expectations from the ensemble tests and a systematic error of 0.2% (2.7%) is assigned 
in the B + — > <f)K +/ ~f (B° — > <pK 0/ -f) mode. The largest contribution comes from the non- 
resonant yield (8.0%). The total systematic uncertainty assigned to the estimated yield 
is 8.2% (8.8%). We also assign a systematic error of 3.3% (4.6%) due to uncertainty on 
charged track efficiency, 1.4% due to particle identification, 2.4% due to photon detection 
efficiency, 1.4% due to uncertainty in the number of BB pairs in B + — > cf)K +r j (B° — > (pK '~f) 
mode. Furthermore, we assign a systematic error of 4.6% in the neutral mode due to Kg 
reconstruction. The statistical uncertainty on the MC efficiency after reweighting is 0.9% 
(1.3%). The uncertainties due to daughter branching fractions account for a systematic 
contribution of 1.2%. We add each contribution above in quadrature to obtain the total 
systematic uncertainty of 9.9% (11.9%). 
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In summary, we report the first observation of radiative B° — > (j)K 0, y decays in Belle 
using a data sample of 772 x 10 6 BB pairs. The observed signal yield is 35 ± 8 with a 
significance of 5.4 a including systematic uncertainties. The measured branching fraction 
is B(B° 0if°7) = (2.66 ± 0.60 ± 0.32) x 10~ 6 . We also precisely measure B(B + -> 
(j)K + -i) = (2.34 ± 0.29 ± 0.23) x 10~ 6 with a significance of 9.6 a. The signal events are 
mostly concentrated at low 4>K mass, which is similar to a two-body radiative decay. The 
neutral mode has enough statistics to measure time- dependent CP asymmetry in order to 
search for new physics from right-handed currents in radiative B decays. 
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